Multiperiodic magnetic structures in Hubbard superlattices 
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We consider fermions in one-dimensional superlattices (SL's), modeled by site-dependent Hubbard- 
U couplings arranged in a repeated pattern of repulsive (i.e., U > 0) and free (U = 0) sites. Density 
Matrix Renormalization Group (DMRG) diagonalization of finite systems is used to calculate the 
local moment and the magnetic structure factor in the ground state. We have found four regimes 
for magnetic behavior: uniform local moments forming a spin-density wave (SDW), 'floppy' local 
moments with short-ranged correlations, local moments on repulsive sites forming long-period SDW's 
superimposed with short-ranged correlations, and local moments on repulsive sites solely with long- 
period SDW's; the boundaries between these regimes depend on the range of electronic densities 
p and on the SL aspect ratio. Above a critical electronic density, p-fi, the SDW period oscillates 
both with p and with the spacer thickness. The former oscillation allows one to reproduce all SDW 
wave vectors within a small range of electronic densities, unlike the homogeneous system. The 
latter oscillation is related to the exchange oscillation observed in magnetic multilayers. A crossover 
between regimes of 'thin' to 'thick' layers has also been observed. 



I. INTRODUCTION 

Magnetic multilayers have been the subject of intense 
study over the last decade. The technologically impor- 
tant giant magnetoresistance (GMR) is one of the most 
interesting aspects of these compounds. Another aspect 
which has brought attention to multilayers is the oscil- 
lation of the exchange coupling between magnetic lay- 
ers as the spacer layer thickness is varied. While os- 
cillations with single periods have been well understood 
for some rljWie, multiperiodicity has been theoretically 
predicted pel and indeed observed, in trilayer materi- 
als. Fe/Cr/Fe samples grown by sputtering or molecular- 
beam epitaxy (MBE) display two periods of oscillation of 
the exchange coupling: a so-called long-period, of about 
10 to 12 monolayers thick, is superimposed to arshort- 
period component of about two monolayers thick.cl This 
superposition of short- and long-period components has 
also been observed in other MBE-fabricated trilayex ma- 
terials such-as Fe/Mn/Fe,t] Fe/Au/Feja Fe/Mo/Feja and 
Co/Cu/Co.0 Short-period oscillations, however, disarm 
pear if interface quality is not carefully maintained.Htm 
Recent experiments!!^ in Fe/Cr/Fe show that areas of 
constant Cr thickness, with diameter larger than 3-4 nm 
on the interface, are necessary for the development of 
short-period oscillations. It is therefore believed that 
multiperiodicity has not yet been observed in multilayers 
due to interface roughness. 

From the theoretical point of view, both the Quan- 
tum well theoryB and the so-called RKKY theoryB can 
account for many features related to the oscillations of 
the exchange coupling. For instance, a direct relation 
between the periods of oscillation and Fermi -surface ex- 
trema of bulk spacers has been established SMB How- 



ever, since the notion of a Fermi surface is not widely 
applicable to strongly correlated systems, a deeper un- 
derstanding of multiperiodicity is clearly in order, and 
microscopic models should provide useful insights. 

With this in mind, here we investigate the mag- 
netic properties of a one-dimensional superlattice (SL) 
modelijii-a in which electronic correlations are incorpo- 
rated and treated nonperturbatively. The model consists 
of a periodic arrangement of Ljj sites ( "layers" ) in which 
the on-site coupling is repulsive, followed by Lq free (i.e., 
U = 0) sites. The role played by relative layer thicknesses 
on the magnetic and conducting properties of these sys- 
tems is then probed by varying Lq and Ljj. 

The SL structure gives rise to several remarkable 
features,!! 2 ] in marked contrast with the otherwise homo- 
geneous system: Local moment maxima can be trans- 
ferred from repulsive to free sites, and the range of pa- 
rameters in which this .occurs has been expressed in terms 
of a 'phase diagram'. lij In addition, spin-density- wave 
(SDW) quasi-order can be wiped out as a result of frus- 
tration, and the SL structure also induces a shift in the 
density pi at which a Mott-Hubbard insulating phase sets 
in.L3 Furtherj-by examining the Luttinger liquid version 
of the model,E3 one finds that these superlattices provide 
the means to realize gapless insulating phasesS3 

Preppus studies of the discrete version of the 
modelEaiiJ resorted to Lanczos diagonalization, which 
sets limits on the system sizes used; for instance, a 24- 
site lattice size could only be considered for the low- and 
high-density regimes (p = 1/6 and p — 11/6). Nonethe- 
less, one was still able to probe the period of exchange os- 
cillations for these special densities through the analysis 
of the magnetic structure factor: the peak position dis^ 
played oscillatory behavior with the spacer thickness. E£l 
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Here we use the Density Matrix Renormalization Group 
(DMRG) technique^-] to study superlattices longer than 
the ones available through the Lanczos method. With 
the aid of the magnetic structure factor, we have been 
able to probe the periodicity of the superlattice over a 
wider range of layer thicknesses and densities. As we 
will see, this has led to significant improvements on the 
phase diagram previously reported,Ej with the addition 
of information relative to the regions in which one- and 
two-period-oscillations are found; as it turned out, these 
regions are closely related to the behavior of the local 
moment. We have also been able to observe a crossover 
between the regimes of thin and thick layers; in the latter 
regime, the 'aspect ratio' £ = Ljj / Lq is the only relevant 
geometric parameter, whereas the magnetic behavior in 
the former regime depends on Ljj and Lq separately. 

The layout of the paper is as follows: In Section y we 
introduce the superlattice model and comment on the 
calculational procedure. Section III focus on the local 



moment and how it changes with density and layer thick- 
ness. The magnetic structure factor and the periodicity 
of the superlattices are discussed in Section IV, and Sec- 
tion |v| summarizes our findings. 
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FIG. 1. Local moment as a function of the site (i) for the 
SL with Lu = 1, L = 2, U = 4, N a = 48, for p = 0.25 
(squares), p — 0.667 (triangles) and half-filling (circles). The 
local moment profile changes qualitatively as the density in- 



II. MODEL AND CALCULATIONAL 
PROCEDURE 

We define the Hamiltonian as 
U = -t ^2 ( c L c i+ia + H.c.J + ^ Ui numx (1) 

where, in standard notation, i runs over the sites of a 
one-dimensional lattice, c ia (ci a ) creates (annihilates) 
a fermion at site i in the spin state a = | or J,, and 
Ui = riii + n iii with n icr — c i(y Ci a \ the on-site Coulomb 
repulsion is taken to be site-dependent: Ui = U > 0, for 
sites within the repulsive layers, and Ui — otherwise. 

We consider the Hamiltonian (|l|) on lattices with N s 
sites and N e electrons, and open boundary conditions 
are used. The appropriate Finite Size Scaling (FSS) pa- 
rameter, however, is the number of periodic cells, N c = 
N s /Nb, for a basis with Nf, — Ljj + Lq sites. The ground 
state wave function and energy are obtained by numeri- 
cal diagonalization using the finite-system-density-matrix 
renormalization group (DMRG) methodO We used lat- 
tice sizes up to 150 sites, and truncation errors in the 
DMRG procedure were kept around 10~ 5 or smaller. We 
have performed a systematic study of the magnetic prop- 
erties for different values of the Coulomb repulsion U, 
different occupation p = N e /N s and different configura- 
tions {Ui}. Not all configurations {Ui} fit into all sizes 
and occupations but, since DMRG allows us to study 
a wide range of lattice sizes, we were able to establish 
overall trends. 



III. LOCAL MOMENT PROFILE 

The local moment at site i is defined as {Sf) = | {{n^ — 
nil) 2 ), an d is a measure of both the magnetism and the 
degree of itinerancy of the system. Figure |l| shows the 
local moment profile for the SL with Ljj = 1, Lq = 2, 
U = 4, N s = 48, and for three different densities; effects 
of system size on the local moment are negligible. For 
small densities, such as for p = 0.25, one identifies small- 
amplitude oscillations in the local moment profile; their 
period (27r/2/ep, with Ikp = np) is determined not by 
the underlying SL structure, but by the Friedcl oscilla- 
tions infihe charge density of the otherwise homogeneous 
system.tl 

As the density is increased, the SL structure dominates 
over the Friedel oscillation as evidenced by the data: for 
p = 0.667 the maxima lie on the free sites and the mod- 
ulation of the profile perfectly matches that of the SL. 
For large enough densities the maxima migrate to the re- 
pulsive sites, as shown by the data for half filling. One 
should also note that even at the maxima (Sf) is consid- 
erably reduced from its value at the completely localized 
limit (U = oo), namely (Sf) =3/4; the itinerant behav- 
ior in these cases is therefore evident. 

The above example illustrates the existence of three 
regions, characterized by different local moment pro- 
files: homogeneous (or Friedel- like), free-site peaked, and 
repulsive-site peaked. In order to locate the boundaries 
between these regimes it is useful to determine how the 
local moment at repulsive and free sites separately change 
with the density. In addition, we define a bias of the local 
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FIG. 2. Bias (triangles) and local moment at repulsive 
(squares) and free (circles) sites, as functions of density, for 
U = 4 and for superlattices with (a) Lu = Lq = 1, N s = 24; 
(b) Lu = 2, L = 1, N s = 48; and (c) Lu = 3, L = 1, 
N s = 64. 



FIG. 3. Bias (triangles) and local moment at repulsive 
(squares) and free (circles) sites, as functions of density, for 
U — 4 and for superlattices with (a) Lu — Lo — 1, N a = 24; 
(b) Lu = 1, Lo = 2, N s = 48; and (c) Lu = 1, L = 3, 
N s = 64. 



moment maxima as 
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and also study its dependence with the density. 

Figure || shows the local moment [both at repulsive 
(iSf,)) and free sites ((-S^))] and the bias as functions of 
the density, for U — 4. In the case of Fig. 0, Lu > £o, 
with all SL configurations having Lo = 1, and Ljj = 1 
(JV S = 24), 2 (JV S = 48) and 3 (N s = 64). In order to 
reduce the effects of open boundary conditions we have 
averaged over the 6 innermost cells. As the density is 
increased from the completely empty system, we see that 
for densities smaller than po, given by 



Po = 



Ln 



(3) 



the local moment increases, and is the same on both sub- 
lattices (hence 5 — 0). For the SL's with Lq = 1 and 
Lu = 1, 2 and 3 one has po = 0.5, 0.33 and 0.25, re- 
spectively, which are indicated by arrows in Fig.|| This 
density corresponds to having one electron on each cell, 
so that in the case Ljj > Lq electrons have equal proba- 
bility of being either on a free or on a repulsive site for 
P < Po- 

From Fig. || one sees that there is a range of densities 
above po, in which the local moment grows slower on the 
repulsive sites than on the free ones, since added electrons 
will preferentially occupy the free sites; hence a negative 
bias develops within this range. By the same token, (Sq) 
will reach its maximum value at densities smaller than 



those at which displays its maximum; for complete- 
ness, recall that the maximum value of the local moment 
on a homogeneous free lattice is 0.375, occurring at half 
filling. 

In strong coupling, the free sites saturate at the density 
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L + L L 



(4) 



which corresponds to having two electrons on each free 
site, while the repulsive site is empty. Nonetheless, even 
for moderate couplings, this density is special. Indeed, 
from Figs. ||(b) and |](c) one can see that for Lu > Lo 
the bias reaches its minimum value exactly at Also, 
the local moment at repulsive sites shows a bump at p^ , 
indicating the beginning of a steady occupation of repul- 
sive sites. 

Increasing the density even further, one sees that (5^) 
reaches its maximum at pi, defined as 



Pi 



2L 



Lq + Lu 



(5) 



which, in strong coupling, corresponds to having two elec- 
trons on each free site and one on each repulsive site; the 
maximum of {S^,), at exactly this density, is indicative 
of the SL being in a Mott-Hubbard insulating state.Ed In 
the region between and pi, the repulsive layer is pref- 
erentially filled as the overall density is increased, causing 
a steep rise (drop) in the local moment at the repulsive 
(free) sites. As Fig. |^ shows, the consequence is a steady 
increase of the bias in this interval. 
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FIG. 4. Magnetic structure factor [Eq. (g)] for a SL with 
Lu = 2, L = 1, N s = 60, N e = 52 (hence p = 0.87), and for 
U = 4 and 16. 



FIG. 5. Magnetic structure factor [Eq. for a SL with 
Lu = 2, Lo = 1, U = 8, p — 1.2, and for different system 
sizes: N s — 30 (squares), N 3 — 60 (circles), N s — 120 (up 
triangles), and N s = 150 (down triangles). 



For densities larger than pi the free sites are almost 
completely filled, which is apparent by a considerable 
decrease in the magnitude of the derivative of (Sq) with 
respect to the density. Fcrmions will then start to dou- 
ble occupy the repulsive sites, thus causing a reduction 
in (S*>. * 

The equivalent of Fig. |2| for the case Ljj < Lq is shown 
in Fig. [| While the overall behavior is the same, a few 
differences are worth mentioning. The first one is the 
behavior of the bias in the range p < p : While the bias 
vanishes for Lu > Lq, for Lu < Lq it is negative, though 
of small magnitude. This is due to the fact that in this 
range of densities, and within each cell, the electrons have 
more free sites at their disposal to resonate than repulsive 
ones; this excess of free sites within each cell also explains 
why the bias still decreases for densities above po. 

Second, for Lu < La, (Sf}) is boosted whenever p = 
2mpo, with m — 1, 2, . . . , Lq + Lu — 1; this can be at- 
tributed to the fact that the double occupancy of the 
repulsive sites is least likely whenever there are an even 
number of electrons per cell. Note also that the first 
bump, at 2po, coincides with the minimum of the bias. 

And, third, while for Lu > Lq the bias changes sign 
for < p < pi, when Lu < Lq this occurs for po < 
p < p^i; the actual location of the density at which 5 = 
depends on the SL configuration, as well as on the 
Coulomb repulsion U. 



IV. MAGNETIC STRUCTURE FACTOR AND 
EFFECTIVE DENSITIES 

Let us now turn to the magnetic structure factor, which 
is defined as 

5(g) = ^^e^-^(Si-Sj). (6) 

As q is related to the repeating units, S(q) probes the 
relative arrangement between different cells. It is im- 
portant to have in mind that the homogeneous system 
displays a single peak in the magnetic structure factor at 
qmax = 2fc Fr-y Trp, for p < 1 , or q max = 2k F = 7r(2 - p), 
for p > l;El the lattice spacing is taken to be unity 
throughout this paper. 

Figure ^ shows S(q) for a SL with Lu = 2, Lq = I, 
p = 0.87, and for two values of U , namely U = 4 and 
U = 16. Two peaks in the magnetic structure factor are 
clearly seen in this case: one at q = it, and another at 
q = 3tt/5. While the former is not affected by an increase 
in U, the latter grows with U, though without changing 
its position. Actually, for sufficiently large U the peak at 
q 7^ 7T even becomes more pronounced than the one at 
q = tt. Further data show that this happens for a range 
of values of Lq, Lu, and p, as discussed below. 

The presence of two peaks (at, say, q max and q' max , 
with q' max < q m ax) in the structure factor is associ- 
ated with a tendency of the system to order (strictly 
speaking, to quasi-oidev, in one dimension) in a mag- 
netic arrangement dominated by the corresponding pe- 
riods, A = 2-K/q max and A' = 2ir/q max . As we will see 



4 



below, the long period oscillates with the spacer thick- 
ness, a behavior reminiscent of the exchange oscillation 
observed in magnetic trilayers. 

These two peaks also differ in the way they depend on 
the system size. Figure |5] shows S(q) for the SL with 
Lu = 2, L = 1, U = 8, p = 1.2 ( Pn < 1.2 < pi), and 
for four different lattice sizes, ranging from N s = 30 to 
N s = 150. From Fig. | we see that the inflection already 
present for N s = 30 at q = 2ir/5 sharpens as N s increases, 
and that there is no change in the position of the peak. 
We have checked that a similar slow, but steady, growth 
of the peak height with N s occurs for the homogeneous 
Hubbard model away from half filling. These features 
have been observed for other SL configurations and den- 
sities, which therefore indicate that whenever a peak is 
found at q ^ 7r, it is robust. On the other hand, the peak 
at q = 7r shows a much weaker size-dependence, so that it 
should be associated with strong, although short-ranged, 
correlations; this point is illustrated below. 

We can then turn to a systematic study of the number 
of peaks and their positions, by analysing the evolution 
of the structure factor as the density of electrons is in- 
creased. As discussed in Sec. Ill, for p < po the local 
moments are small, and either their maxima are on the 
free layers, or they are evenly distributed throughout the 
lattice, depending on whether Ljj < Lq or Ljj > L$, re- 
spectively. A small or zero bias signals that the SL struc- 
ture is not very relevant in this situation. Indeed, the 
spatial decay of the spin-spin correlation function (not 
shown here) in the case of a superlattice with a small 
bias can hardly be distinguished from that of the corre- 
sponding homogeneous system; as a result, the magnetic 
structure factor displays a single peak. In addition, this 
single peak displays a size- and a [/-dependence similar 
to those for the homogeneous system. 

In order to relate the peak position with some density, 
one can think of a free (homogeneous) lattice in which 
the sites are grouped in cells mimicking the SL structure 
under consideration; it then follows that a meaningful 
quantity is the cell density of electrons, 
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p{L v + L ), 



(7) 



where p is the overall density. For the interacting SL, we 
have found that the peak position is given by the same 
expression as for the homogeneous case, but with p ce ii 
replacing p; that is, 



TTPcelh f°r P < Po- 



(8) 



Thus, the peak position grows linearly with p up to p — 
Po (at which density p cc u = 1), when it reaches q m ax = R] 
see Figs. || and |?]. 

The single peak regime persists for po < P < Pu> an d> 
as shown in Figs. |(] and Q now the peak is always at 
q = it. The single peaks in this region show a very weak 
dependence on the system size, which is reflected in the 
spatial decay of the correlation functions, (SqSj). As 
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FIG. 6. Maxima position, q m ax, of the magnetic structure 
factor as a function of density for U = 4, and (a) Lu = 1 
L = 1 N s = 24, (b) Lv = 1 L = 2 N s = 48, and (c) Lu = 1 
Lo = 3 N s = 64. The dashed lines indicate the presence of 
another peak in S(q) (see text). 



illustrated in Fig. ^(a) for Lu = Lo = 1, U = 4, and 
p = 0.75, correlations with origin in either of the sublat- 
tices barely survive at large distances; this should be con- 
trasted with the case displayed in Fig. |8|(b), for p = 1.75 
(see below), in which correlations in one of the supcrlat- 
tices are 'long' ranged. 

At pfj_, and in strong coupling, the free layers are com- 
pletely filled while the repulsive layers are empty. But 
as the density is increased beyond pn, a second peak 
emerges, as indicated by the dotted lines in Figs. and 
0. This second peak results from the robust moments 
located on the repulsive sites. Indeed, if one defines an 
effective electronic density on the repulsive layers as 



Pett = p(Lo + Lu) - 2L , 



(9) 



where p is the overall density, the long-periods are located 
at 



( 1„ 



(10) 



With this definition, it also becomes clear that for p = 
Pll there is no net moment at the repulsive layers, since 
Pcff = 0. 

This two peak structure is present until one reaches 
pu, defined by 



pu = 2- po, 



(11) 



which corresponds to a single hole per cell. One should 
also have in mind that the overall magnetic arrangement 
is determined by the long-period (characterized by q' max ), 
since, as discussed above, this is the one increasing with 
system size. 
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FIG. 7. Maxima position, q m ax, of the magnetic structure 
factor as a function of density for U — 4, and Ljj = 1 Lo = 1 
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Lo = 1 iV s = 64 (c). The dashed lines indicate the presence 
of another peak in 5(g) (see text). 
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FIG. 8. Spatial decay of correlations, for a SL with 
Lv = L = 1 and U = 4 for N s = 48 sites: in (a) p = 0.75, 
and in (b) p = 1.75. Circles and squares correspond to the 
origin being taken on free and repulsive sites, respectively. 



An interesting difference between the cases depicted in 
Figs. | and @ is the fact that in the former pjj = pi, while 
in the latter pjj > pi, and q' max is able to go through at 
least one complete oscillation before p reaches pjj. In 
this case, the situation q' max = does not indicate any 
tendency towards a ferromagnetic arrangement, but is to 
be associated-arith frustration of the corresponding long- 
period SDW.El Indeed, when L v = 2 and L = 1 the 
Mott-Hubbard insulator at pj is frustrated, since two 
spins on each repulsive layer form local singlets. Sin- 
glets on different repulsive layers, in turn, do not cou- 
ple with each other, though short ranged correlations are 
still present; see Fig. |^(b). The frustration at half fill- 
ing for Ljj — 3 and Lq = 1 can be understood by a 
similar strong coupling analysis: of the four electrons on 
each cell, two occupy the free site and the remaining two 
resonate between three sites, but always forming a sin- 
glet. Figure 0(c) shows that further addition of electrons 
renders these singlets unfavorable, and the system again 
displays a SDW. At p = 3/2, one reenters a frustrated 
state, again as a result of having an even number of elec- 
trons on the repulsive layer. Therefore, we can relate 
the reentering frustrated configurations to the formation 
of singlets on the repulsive layer, which occurs whenever 
there is an even number of electrons per cell; that is, 
whenever the density goes through an even multiple of 
2po- 

And, finally, above pjj all SL's return to a single peak 
regime: S(q) has a maximum at 7r(2— p ce u) [= 7r(2 — p e g), 
since p c g — p cc \\ — 2L ]. The correlations in this regime 
are quasi-long-ranged, since, as shown in Fig. ||(b), the 
correlation function with origin at a repulsive site is 
slowly decaying. 

The above analyses of the magnetic structure factor 
and of the local moment profile can be extended to sev- 
eral other SL configurations, and the outcome can be 
best summarized by a diagram in the parameter space 
(p, Lq, Ljj), showing the presence of four different regions 
(or phases). Cross sections of the full three-dimensional 
phase diagram are presented in Fig. ^|(a) for Ljj = 1 
and Fig. |](b) f° r Ljj = 4. In the low density region 
(A), located between p = and p = p , the system 
behaves roughly as if it were homogeneous. The local 
moments are small and their maxima are located prefer- 
entially on the free layers. The SDW is dominated by a 
single density-dependent wavevector, q max = npceii- 

At somewhat larger densities, po < p < p-^, lies a 
region (B), in which the positions of the maxima in the 
local moment profile depend on the repulsion U, on Lq, 
and Ljj- Presumably as a result of this 'floppy' character 
of the local moments, spin correlations in this region are 
strongly antiferromagnetic, but short-ranged. 

As the density is further increased, one enters the 
double-period region (Ci), which lies between p = p-jq 
and p = pjj. In this region the local moment on the 
repulsive layer suffers successive boosts, and one finds 
SDW's with a 'long' period A' = 2/p e s; the latter are ac- 
companied by strong short ranged correlations, of period 
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FIG. 9. Regions in the parameter space (p,Lo) for (a) 
Lu — 1 and (b) Lu = 4: A - weak moments formed pref- 
erentially on free layers and one single SDW period; B - lo- 
cal moment maxima depend on U and on Lq/Lu and spin 
correlations are predominantly antiferromagnetic, but short 
ranged; CI - local moment maxima on the repulsive layers 
and the SDW's are dominated by two periods; C2 - local 
moment maxima on the repulsive layers and the SDW's are 
dominated by a single period. See text for details. Full line 
corresponds to po, dotted line to pfi and dashed line to pu- 



A = 2. 

And, finally, there is a high density region (C2), with 
densities ranging from p = pu to p = 2. At pu the 
local moment bias is maximum (see Figs. |^ and ||), so 
it decreases as one increases the density. Nonetheless, 
one still has SDW's, now with a single period given by 
A = 2/(2 — Pccii)- By comparing the two cases depicted 
in Fig. y, one sees that a growth of the repulsive layer 
increases the two-peaked region, at the expense of all 
others. 

The full three-dimensional phase diagram is shown in 
Fig. [Io| The densities po, p^i, and pu define surfaces in 
the parameter space (p, Lq, Lu) which act as boundaries 
between the four regions discussed above. The po-surface 
flattens considerably for thick layers, and if one imagines 
an £ = 1 line on the horizontal plane of the figure, we 
see that the homogeneous-like region is only important 
for moderately thin layers. The pjj-surface is the same 
whether the layers are short or long, since it depends on 




FIG. 10. Regions in the parameter space (p, Lo,Ln): the 
lower surface corresponds to po, the middle one to pu, and 
the upper one to pu- The dotted line is the intersection of 
the p|j.-surface with the plane £ = 1 (see text). 



Lu and L only through the combination Lu / L = I; 
to illustrate this, the intersection of the pti-surface with 
the plane 1 = 1, shown as a dotted line in Fig. [n], yields 
PU = 1 for all Lu = Lq. The topmost surface (pu) 
also displays a similar crossover between thin and thick 
regimes: for thick lattices pu — > 2. 

We are now in a position to discuss the oscillation 
in Qmax with the spacer (free layer) thickness, for a 
fixed electron density; as mentioned before, these are re- 
lated to the oscillation of the exchange coupling between 
magnetic layers. When p < po, one has q m ax < it, 
so that there is no oscillatory behavior in q max . For 
Po < P < PTJ.1 the peak is always at q max — n, so 
that again no oscillation is found. Above the long- 
period maxima in the magnetic structure factor are lo- 
cated at q m ax(Lo) — 7rp e ff(Lo), where we have empha- 
sized the dependence with Lq through p e g. We can 
then calculate the period of oscillation, ALq, by setting 
(Lo) (L + AL ), mod(27r). For p n < p < 1 

we get 



AL =(1-^ 

7T 



(12) 



where 2k p — np since p < 1; for 1 < p < pu , we similarly 
find 



AL 



(13) 



where now 2kp = ir(2 — p). Note that this result is not 
valid for Lu — 1, since, according to Fig. 0, once q max 
vanishes, it does not grow as Lq increases. For p > pu, 
on the other hand, Eq. ([l3]) is applicable. 

We have then established that (i) p^ acts as a criti- 
cal density for the appearance of 'exchange oscillations', 
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and that (ii) our previous results^ for ALq, obtained 
in the high density region, are valid quite generally for 
p > /9|j_. Further, Eqs. ( [l2"| ) and ( |l3| ) reproduce previous 
findings, within the Hartree-Fock approximation, for the 
periods of oscillation of the exchange coupling in mag- 
netic multilayers. ErEil Also, the experimentally observed 
short period of two monolayers reported in Ref. || cor- 
responds, in our framework, to the A = 2 correlations. 
Thus, electronic correlations do not modify the quantum 
interference effects determining the periods of oscillation 
from the extrema of the Fermi surface of the spacer ma- 
terial. 



V. CONCLUSIONS 

We have investigated one-dimensional Hubbard super- 
lattices consisting of periodic arrangements of free and 
repulsive layers, by means of Density Matrix Renormal- 
ization Group. By considering a much wider range of lat- 
tice sizes and densities than in previous studies, we have 
refined in several aspects our earlier predictions for the 
magnetic behavior. There are now four distinct regimes, 
depending on the range of electronic densities. For less 
than one electron per periodic cell, the local moment pro- 
file is approximately uniform, and spin-density waves are 
dominated by a single density-dependent wave vector. 
When the density lies between those corresponding to 
one electron per cell and to a fully occupied free sublat- 
tice (with empty repulsive sites), maxima in the local mo- 
ment profile develop, which can be either on free sites or 
on the repulsive sites, depending on the SL configuration, 
on the density, and on U; also, spin correlations become 
short ranged, but dominated by a tendency of neighbor- 
ing cells to align antiparallel. For densities larger than 
two electrons per free site one has a two-period magnetic 
structure. There is a long period SDW, in which the 
wave vector oscillates as a function of the electronic den- 
sity. An immediate consequence is that SDWs with all 
possible wave vectors are generated within an interval of 
densities of 2po; this should be compared with the homo- 
geneous system, for which one needs to vary between an 
empty lattice and a half filled one in order to generate 
all possible wave vectors. These long-period SDWs are 
superimposed with short ranged correlations with q = ir, 
which disappear for densities above one hole per periodic 
cell. 

We have also extended to a broader range of densities 
our earlier prediction that the wave vectors for the SDWs 
oscillate as the free layer length is varied, with a period 
determined solely by the electronic density (through the 
Fermi wave vector). In the context of magnetic multi- 
layers, our results for the period of oscillations exactly 
reproduce the relation between Fermi surface extrema 
with exchange coupling oscillation; in addition, the two- 
monolayer period observed experimentally corresponds 
in our model to the short period at q = it. 



And, finally, we have been able to observe a crossover 
between the regimes of thin and thick layers; in the latter, 
the 'aspect ratio' £ = Ljj/Lq is the only relevant geomet- 
ric parameter, whereas in the former regime the magnetic 
behavior depends on Ljj and Lq separately. For instance, 
when any of the layers are thin - less than about 6 sites 
long -, the SL structure is not felt at low densities, and it 
behaves as if it were homogeneous; for thick layers, this 
quasi-homogeneous behavior is only noticeable at very 
low densities. Similarly, the region of singly peaked cor- 
relations at high densities gets smaller as the layers get 
thicker. 

As a final comment, one should expect that the ap- 
plicability of the one-dimensional model treated here 
is very close to being extended beyond the realm of 
higher dimensional superlattices. Indeed, fabrication of 
nanowire superlattices has been recently reported.E2l Al- 
though these superlattices were made up of semiconduct- 
ing materials, the prospects of growing metallic and/or 
magnetic nano superlattices are promising. In this case, 
our results indicate that a careful control of the doping 
level leads to a wide variety of distinct magnetic behav- 
iors in the same material. Another possible realization of 
our model would be to a (as yet hypothetical) superlat- 
tice made up of single- walled metallic carbon nanotubes, 
since these have been successfully described in terms of 
a Luttinger liquid; see, e.g., Ref. 16 for a partial list of 
references. 
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